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Amplitude damping changes entangled pure states into usually less-entangled mixed states. We show, how- 
ever, that even local amplitude damping of one or two qubits can result in mixed states more entangled than 
pure states if one compares the relative entropy of entanglement (REE) for a given degree of the Bell-CHSH 
inequality violation (referred to as nonlocality). By applying Monte-Carlo simulations, we find the maximally- 
entangled mixed states and show that they are likely to be optimal by checking the Karush-Kuhn-Tucker condi- 
tions, which generalize the method of Lagrange multipliers for this nonlinear optimization problem. We show 
that the REE for mixed states can exceed that of pure states if the nonlocality is in the range (0,0.82) and the 
maximal difference between these REEs is 0.4. A former comparison [Phys. Rev. A 78, 052308 (2008)] of 
the REE for a given negativity showed analogous property but the corresponding maximal difference in the 
REEs is one-order smaller (i.e., 0.039) and the negativity range is (0,0.53) only. For appropriate comparison, 
we normalized the nonlocality measure to be equal to the standard entanglement measures, including the nega- 
tivity, for arbitrary two-qubit pure states. We also analyze the influence of the phase-damping channels on the 
entanglement of the initially pure states. We show that the minimum of the REE for a given nonlocality can be 
achieved by these two channels, contrary to the amplitude damping channel. 



PACS numbers: 03.65.Ud, 03.67.Mn 

I. INTRODUCTION 

Quantum entanglement and nonlocality yj, |2[] are among 
the central concepts in quantum theory providing powerful 
resources for modern quantum-information processing [jH. 
Nevertheless, despite of much progress (as reviewed in, e.g., 
Ref. J4[]), our understanding of the relation between entangle- 
ment and nonlocality is very incomplete. 

Quantum nonlocality is often discussed in the context of 
Bell's theorem, which can be, arguably, considered „the most 
profound discovery of science" [5]. The Bell inequality vio- 
lation (BIV) for a given state implies that there is no physi- 
cal theory of local-hidden variables which can reproduce the 
predictions for the state IUJ . Although the BIV does not nec- 
essary imply nonlocality and vice verse' a (see, e.g., Ref. Ja]), 
for convenience, we use the terms BIV and nonlocality inter- 
changeably. 

Because of the prominent role of quantum entanglement, 
much effort has been devoted to investigating measures of en- 
tanglement that can be used for its quantification. Currently, a 
variety of formal and operational entanglement measures can 
be applied |4j] depending on the physical and information con- 
texts. For example, if one is interested in distinguishing entan- 
gled from separable states, a convenient choice is the relative 
entropy of entanglement (REE) [7]. However, if one considers 
the cost of entanglement under operations preserving the pos- 
itivity of partial transpose (PPT) [8], the appropriate choice of 
entanglement measure is the (logarithmic) negativity Ir9i [Toll . 
The concurrence 111 ill is yet another popular measure of entan- 
glement, which quantifies the entanglement of formation Il2ll . 
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In the case of two-qubit pure states some of the entangle- 
ment measures are equivalent (e.g., the negativity, concur- 
rence and nonlocality). However, the situation is quite dif- 
ferent for mixed states, where the states of the same value 
of a given entanglement measure can correspond to different 
values of entanglement quantified by other measures. Thus, 
it can happen that entangled states ordered according to one 
entanglement measure can be ordered differently by another 
measure lfl3l - tl6ll . 

The most interesting entangled states are, arguably, the 
boundary states, which are extremal in one measure for a 
given value of another measure. They can serve as "way- 
marks" for entangled states providing means for their system- 
atic classification. Moreover, it turns out that some of the most 
entangled states according to the REE for a given value of an- 
other entanglement measure as, e.g., the negativity 01711 . are 
mixed states. 

Some comparative studies of the nonlocality in relation to 
entanglement measures in a dynamical context of two qubits 
can be found in, e.g., Refs. 03 SH- 

In this paper, we find maximally entangled mixed states 
(MEMS) corresponding to the largest REE for a given degree 
of nonlocality using the Horodecki measure ll2lll22ll . We also 
show that these MEMS are more entangled than pure states 
for a wide range of this nonlocality degree. As an indicator 
of extremality of the found MEMS we use the Karush-Kuhn- 
Tucker (KKT) conditions in a refined method of Lagrange 
multipliers. 

The paper is organized as follows. In Sec. [II] we intro- 
duce some basic definitions of the nonlocality and entangle- 
ment measures used in our paper. In Sec. [Hi] we demonstrate 
analytically that there are mixed states more entangled than 
pure states when analyzing the REE for a fixed nonlocality 
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(and negativity). In Sec. IIV1 we describe how to physically 
obtain an important class of entangled states, which, as we 
show in Sec. [V] are the boundary states if quantified by the 
REE for given values of the nonlocality, negativity, and con- 
currence. The most important result in this paper is finding 
that the states exhibiting the highest and lowest REEs for a 
given value of the nonlocality and demonstrating that mixed 
states can be more entangled than pure states if the fixed non- 
locality is less then 0.8169. 

II. DEFINITIONS 

Hereafter, we study general two-qubit density matrices p, 
which can be expressed in the standard Bloch representation 

as 

1 3 
p = -(I®I + x-a®I + I®y-a + ^2 T mn a n <E)<T m ), (1) 

n.m— 1 

where S = [a\, 03] and the correlation matrix = 
Tr[p(cTj <g> CTj)] are given in terms of the three Pauli matri- 
ces, and Xi = Tr[p(<Ti ® I)] (y, = Tr[p(7 ® Oj)]) are the 
elements of the Bloch vector x (y) of the first (second) sub- 
system. Expressing the two-qubit density matrix by Eq. ([T) 
is very convenient since it allows a direct application of an 
effective criterion for the nonlocality. 

A. Nonlocality measure 

The two-qubit Bell inequality in the form derived by 
Clauser, Home, Shimony and Holt (CHSH) |2fl can be for- 
mulated as 

|Tr(pB C HSH)|<2, (2) 

where the Bell-CHSH operator is Schsh is given by 

Schsh = o- <x<g) (b + b') -a + a' ■ a® (b- V) ■ a, (3) 

and its expected value is maximized over real-valued three- 
dimensional unit vectors a, a', b, and b'. According to the 
Horodecki theorem, the maximum expected value of the Bell- 
CHSH operator for a given state p reads as i2~ll I22I1 : 

max Tr (p^chsh) = 2 y/M(p) (4) 

OCHSH 

given in terms of the parameter 

M (p) = max {hj + h k } < 2, (5) 

j<k 

where hj (j — 1,2,3) are the eigenvalues of the real symmet- 
ric matrix U = T T T constructed from the correlation matrix 
T and its transpose T T . Hence, the condition for violating 
the Bell-CHSH inequality is M(p) > 1 JUIH. For con- 
venience, we refer to nonlocality as the violation of the Bell- 
CHSH inequality. 



To quantify a degree of the nonlocality, one can directly use 
M(p) or 2 y/M(p) (see, e.g., Refs. jUB), or more natu- 
rally max [0, M(p) — 1] (see, e.g., Ref. |l8ll ). However, we 
decided to use another function of M(p), denoted by B{p), 
which for two-qubit pure states is equal to the concurrence 
and negativity. This measure can be given as [ 16]: 

B{p) = Vmax[0, M(p)-1]. (6) 

We see that B(p) = if a given state p satisfies the Bell- 
CHSH inequality, given by Eq. @, and B(p) = 1 if the in- 
equality is maximally violated. The value of B > increases 
with M, thus, it can be used to quantify the BIV. We refer to 
this BIV degree as nonlocality measure. 



B. Entanglement measures 

Now, we recall some definitions of a few selected measures 
of entanglement applicable for two-qubit entangled states. 

In our considerations, the most important entanglement 
measure is the REE defined as Er(p) = min (Te x>S(/5||cr), 
which is the relative entropy 5(/o||ct) = Tr (p\o^ p — 
p log 2 a) minimized over a set V of separable states a Il7l l25ll . 
This measure is a quantum counterpart of the Kullback- 
Leibler divergence quantifying the difference between two 
classical probability distributions. Evidently, the REE is lim- 
ited, by definition, to distinguishing a density matrix p from 
the closest separable state (CSS) a only. Note that the REE 
is not a true metric, since it is not symmetric and does not 
fulfill the triangle inequality. However, the REE has a desir- 
able property of a good entanglement measure that, for pure 
states, it reduces to the von Neumann entropy of one of the 
subsystems. 

Unfortunately, as discussed in Refs. 1I26U29I1 . it is very un- 
likely to find an analytical compact formula for the REE of 
a general two-qubit mixed state, which would correspond to 
finding its closest separable state (CSS). Numerical proce- 
dures for calculating the two-qubit REE correspond usually 
to a n op timization problem over 79 or more real parame- 
ters Ill7ll25ll30ll . On the other hand, there is a compact-form 
solution of the inverse problem: If a CSS is known then all 
the entangled states (having the same CSS) can be given ana- 
lytically not only for two qubits ll27il28[ | but even for arbitrary 
multipartite states of any dimensions Bill . 

The second measure studied here is the negativity lfT3[ |32[ 
l33ll . defined as N(p) — max{0, — 2/i in ;„}, where /z m i n = 
mineig( j o r ) and V denotes partial transpose. The negativity 
is related to the logarithmic negativity, log 2 [N(p) + 1], which 
is a measure of the entanglement cost under the PPT opera- 
tions (US]- These two related measures reach unity for the 
Bell states and vanish for separable states, however, for clarity 
of our further presentation, we will use only the negativity. 

The last entanglement measure applied in this paper is 
the concurrence introduced by Wootters 111 ill as C(p) — 
max{0,2A max -X)j A,}, where A| = e\g[p(o- 2 ®a2)p*{p2® 
<T2)]j and A max = maxj Xj. This measure is directly related 
to the entanglement of formation, Ep(p) lfl2ll . However, for 
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the same reason as in the case of the negativity we use the 
concurrence instead of Ep(p). 

III. ANALYTICAL COMPARISON OF ENTANGLEMENT 
AND NONLOCALITY FOR SPECIAL STATES 

In Fig. Q] we presented several curves corresponding, in 
particular, to the Horodecki states, Bell-diagonal states and 
pure states. The REE for these states can be calculated analyt- 
ically so let us first discuss and compare them to demonstrate 
the main point of our paper. 

A. Pure states 

We can simply relate all the above-mentioned entanglement 
and nonlocality measures in a special case of an arbitrary two- 
qubitpure states \ip) = a\00) + b\0l) + c|10) + <2|11) (with 
the normalization condition \a\ 2 + \b\ 2 + \c\ 2 + \d\ 2 — 1) as 

B(\1>)) = C(|V)) = N(\r/j)) = 2\ad - bc\. (7) 

Moreover, in this special case B, N and C are simply related 
to the REE and von Neumann's entropy S as 

W(B) = E R {\i>)) = S(jh), (8) 

where W{B) = h{\[l+ \/l - B 2 ]) is the Wootters func- 
tion 111 ill given in terms of the binary entropy h{x) = 
—xlog 2 x — (1 — x)log 2 (l — x), and pj = Tr 3 _ i p is the 
reduced density matrix of the zth qubit. 

B. Horodecki states 



where B 6 = B(p 6 ) = 0.5856 and Eg >(B 6 ) = Eg '(Be) = 
0.4520 [see Tab. Hand Fig.[3jb)]. This means that mixed states 
can be more entangled than pure states at least for B < B(p§). 
Fig-Q2 c ) shows that this property holds up to B < B(p§) — 
0.8169 but for other mixed states, which is demonstrated in 
El Here, E^ ] as a function of B is explicitly given by 
Eq. (HB forp = ^/(l + B 2 )/2 > l/y/2. 

Analogous comparison of the REEs for a given negativity 
TV for the Horodecki and pure states [see Fig. [TJb)] shows 
that JH: 

^4 H) ( N ) > E P ( N ) for < N < N i i ( 1 3) 
E^p {N) < E { p (N) for JVi < N < 1, 

where ATj ee N( Pl ) = 0.3770 and e£\Ni) = S^ P) (iVi) = 

0.2279 [see Tab. I and Fig. Ha)]). Note that E { ^ } as a 
function of TV can be given explicitly by Eq. ( fTTT i as p = 
V / 2N(1 + N) - N. 

It is convenient to introduce a parameter AEr(X), which 
is the maximal difference in the REE between a given mixed 
state p and some pure state p^ = IV-'XV'I having the same 
value of either the nonlocality (X = B) or negativity (X = 
N), i.e., 

AE R (X') = mzx{E R [p{X)} - E R [p^(X)}} 

= E R [p(X')]-E R [p( p \X% (14) 

where X' is the optimal value of X. For the Horodecki state, 
we observe that this maximal difference for the nonlocality is 
equal to AE R H) (B') = 0.2949, which occurs for B' = 0, 
while for the negativity is only AE^ (N') = 0.0391, which 
is for N' = 0.1540. 



Let us also analyze a mixture of a Bell state, say = 
(|01) + |10))/V2, and vacuum state, i.e., 

P (H) (P) = P (A) (1/2,P) =P|^ + X^ + | + (1-P)!00><00!, (9) 

which is referred to as the Horodecki state. By applying 
Eq. ©, we find that M(p (ii ^) = p 2 + max[p 2 , (1 - 2p) 2 }. 
Thus, the nonlocality B is given by 

S(p (H) ) = Vmax(0,2p 2 - 1). (10) 

On the other hand, the entanglement measures are the fol- 
lowing: the concurrence is C(pW) = p, the negativity is 
N(pW) = -p) 2 +p 2 - (1 - p), and the REE reads as 

E R (pW) = 2h(l - p/2) - h(p) - p 
= (p-2) log 2 (l - p/2) + (1 - p) log 2 (l - p). (11) 

It is seen that the Horodecki state is entangled for < p < 
1, while it is violating the Bell-CHSH inequality only for 
l/y/2 < p < 1. By comparing the REEs for a given non- 
locality for the Horodecki and pure states we find that 

E R H) (B) > E { P{B) for < B < B 6 , (12) 
E ( R n) (B) < E ( R ] {B) for B 6 < B < 1, 



C. Bell-diagonal states 

Finally, let us also analyze the Bell-diagonal states (labeled 
by D), which are defined by 

4 

p^ =£>I/Wil> (15) 

i=l 

where \f3i) are the Bell states and < \j < 1 such that 
J2j = 1 with the largest eigenvalue maxj Xj = (1 + 

N)/2 > 1/2. The nonlocality of the Bell-diagonal state p^ 
is given by [ 16]: 




max{0, 2 max [(A* - A 3 ) 2 + (A fe - A 4 ) 2 ] - 1}, 



(16) 

where subscripts (i,j,k) correspond to the cyclic permuta- 
tions of (1,2,3). By contrast, the negativity and concurrence 
are the same and simply given by N(p^) — C(p^) = N. 
The REE versus the negativity (and, thus, also the concur- 
rence) reads as 

E(p^) = l-h(^) (17) 
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TABLE I. Nonlocality and entanglement measures of the amplitude- 
damped states p n = p' A ' (a,p) (for n = 1, 6), given by Eq. d25t . 
corresponding to the characteristic points marked in Fig. [TJ 



State 


a 


P 


c 


N 


E R 


B 


Pi 


0.0369 


1.0000 


0.3770 


0.3770 


0.2279 


0.3770 


P2 


0.0751 


1.0000 


0.5271 


0.5271 


0.3847 


0.5271 


P3 


0.2198 


0.8536 


0.7070 


0.5756 


0.4039 


0.0000 


Pi 


0.3510 


0.9565 


0.9130 


0.8706 


0.7445 


0.8169 


PS 


0.2116 


1.0000 


0.8169 


0.8169 


0.7445 


0.8169 


P6 


0.0947 


1.0000 


0.5856 


0.5856 


0.4520 


0.5856 



as given in terms of the binary entropy h. If max., Xj < 1/2 
then the state is separable E(p^) = 0. 

It is evident that the nonlocality of the Bell-diagonal states 
depends on all probabilities \, while the entanglement mea- 
sures depend solely on the largest value max^ A^ > 1/2. Nev- 
ertheless, in some special cases of these states, the entangle- 
ment and nonlocality measures can be equal. For example, 
when only two probabilities Ai corresponding to |V' ± ) are 
nonzero, the Bell-diagonal state p^ reduces to 

p P») = p |^+)^+| + (i -p)|Vr><Vr|, (18) 

which is studied in a physical context in Sec. [IV] The nonlo- 
cality for this rank-2 state is simply given by 

B(pW) = \2p-l\ (19) 

implying that 

E( P ^) = l-h(±±^), (20) 

which corresponds to Eq. (fTTt . 

It is seen in Figs. [Tib, c) (see also Sec.[V]for a partial proof) 
that the lower bounds of the REE vs negativity and the REE vs 
nonlocality correspond to the Bell-diagonal state p^°\ which 
satisfy the extremal KKT conditions as we show in Sec. [Vl 
Their physical context is discussed below. 



IV. MANIPULATING ENTANGLEMENT AND 
NONLOCALITY VIA DAMPING CHANNELS 

In Fig. [TJ we have presented the entanglement and non- 
locality measures for a milion randomly-generated two-qubit 
states using the Monte Carlo simulation. It is seen that pure 
states are maximally entangled in terms of the REE for an 
arbitrary concurrence as shown in Fig. [TJ a )- However, pure 
states are not always maximally entangled for the cases shown 
in Figs.[Tfa,b). Namely, the red regions correspond to mixed 
states having the REE higher than that for pure states forgiven 
negativity [Fig. 02b)] and nonlocality [Fig.QJc)]. 

This is a counterintuitive result, so to have a deeper phys- 
ical meaning of the states corresponding to the red regions 
in Fig. QJ let us analyze the loss of entanglement between 




0.2 0.4 0.6 0.8 



Nonlocality 

FIG. 1. (Color online) Monte Carlo simulations of about 10 6 two- 
qubit states p for the relative entropy of entanglement, Er(p), ver- 
sus: (a) concurrence C(p), (b) negativity N(p), and (c) nonlocality 
B(p) corresponding to the Bell-CHSH inequality violation. Red re- 
gions correspond to mixed states more entangled than pure states in 
terms of the En(p). Key: P - pure states, D - the Bell-diagonal 
states (a = 0.5), A' (A") - the MEMS, which are the amplitude- 
damped states optimized for the case b (c), H - the Horodecki states, 
and the coordinates of points p n (where n = 1, 6) are given in 
Tab.|TJ 
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Nonlocality 



FIG. 2. (Color online) Ranges of the relative entropy of entangle- 
ment for a given nonlocality for the two-qubit mixed states gener- 
ated from pure states by the amplitude damping channel, given by 
Eq. d25t (yellow and green areas) and the phase damping channel, 
given by Eq. d32t (blue and green areas). 



two qubits initially in a pure state. In particular, by tak- 
ing one of the Bell states |i/' ± ) = (|01) ± |10))/V2 or 
|0 ± ) = (7(g) eri)|V> ± ) as an initial state, one can achieve any 
degree of entanglement between the two qubits by coupling 
them to their environment(s). Note, however, that by cou- 
pling qubits to separated and uncorrelated environments, we 
can only decrease the entanglement, since any local operation 
cannot increase it. Damping can occur in various ways, but we 
focus on states generated from pure states by three prototype 
damping models (channels), for which the ranges of the REE 
for a given value of the nonlocality is shown in Fig. [2] We 
show in detail that these damping channels provide us with 
the mixed states with extreme values of the REE for a fixed 
value of another entanglement or nonlocality measures. 

In this section we analyze the effect of the prototype damp- 
ing channels on pure states of the form 



|Va) = V^|01) + VT^ a|10), 



(21) 



where < a < 1. When analyzing the entanglement mea- 
sures, any pure two-qubit state, as given above Eq. (0, can 
be equivalently expressed by means of the Schmidt decompo- 
sition, in the form given by Eq. d2"Tl >. Although damping of 
these pure states can lead to inequivalent states, for simplicity, 
we confine our analysis only to the input states \ip a ). 

A convenient description of damping can be given in terms 
of the Kraus operators E% (specified below). Two-side damp- 
ing of the state p- m (a) — \ip a )(ipa\ leads to the following 
output state 

= Y} E Mi)®E 3 {q 2 )]p- m {a) [E\ ( qi )®E}(q 2 )}, 

i,i 

(22) 

where the Kraus operators satisfy the normalization relation 
J2iE}(q)Ei(q) = I. In the special case of the one-side 
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FIG. 3. (Color online) Ranges of the relative entropy of entanglement 
for given values of (a) negativity and (b) nonlocality. It is seen that 
the maximal difference in the REE between the optimal states and 
pure states is equal to (a) A = AEr(N') = 0.0391 at the negativity 
N' = 0.154 and (b) A = AE R {B') = 0.404 at the nonlocality 
B' — 0. Pure states become more entangled than the Horodecki 
states p (H) {p) = p [A \\,p) for (a) N > N(pi) = 0.377 and (b) 
B > B(pe) — 0.586. While pure states become optimal for (a) 
N> N = N(p 2 ) = 0.527 and (b) B > B = B(p 4 ) = B(p 5 ) = 
0.817. The states p n (for n = 1, 6) are specified in Tab.|I] 



damping (say q 2 = 0), Eq. ( 1221 ) reduces to 

p(a, qi ) = ^[^(91) ® I]pin(a)[El( qi ) ® /]. (23) 

i 

A. Amplitude-damping channel 

The Kraus operators for the single-qubit amplitude- 
damping channel (ADC) are the following [[j] : 

£()(%) = |0)(0| + VS|1)(1|, £i(ft) = Vftl°><l|. (24) 

where q t = 1 — p t is the amplitude-damping coefficient. One 
can find that the pure state | t/j a ) after passing through the ADC 
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is changed into the following mixed state 

(a, qi,q 2 )= p\^ ) | + g|00) <00| = (a',p) , 

(25) 

where the effective damping constant reads as q = 1 — p = 
aq 2 + (1 — a)qi and \ip a ') is given by Eq. ( ETT i but for 
q' = ap2/[ap2 + (1 — which can take values from 

to a, where a' = a is achieved for the symmetric two-side 
ADC, i.e., when qi = q 2 . The state given by Eq. ( 1251 ) is some- 
times referred to as the generalized Horodecki state (see, e.g., 
Refs. Hzl |H). Hereafter, we refer to the discussed states 
p( A \a',p) as the amplitude-damped states. The range of 
the REE for a given nonlocality for the ADC state, given by 
Eq. (f25j), is shown in Fig. [2] 

In the special case, when ap 2 = (1 — a)pi, which implies 
a' = 1/2, the ADC reduces to the Horodecki state (see e.g. 
Ref. i e., 

P (H) (P) =P (A Hhp) =P\i> + )0> + \ +9|00)(00|, (26) 

as given by Eq. These states can be obtained from the ini- 
tial Bell state subjected to the symmetric two-side ADC, 
but they can be also obtained by asymmetric two-side ADC 
(or even the one-side ADC) of the initially non-maximally en- 
tangled pure state \tp a ). 

We find that the nonlocality of the ADC state for any a is 
given by 

B(p (A) ) = v/maxjO, max[a;, (1 - 2p) 2 } + x - 1}, (27) 

where x — 4(1 — a')a'p 2 . By comparison, the negativity 
reads as 

N(p^) = ^{1-pf+x - (1 - p), (28) 

and the concurrence is C (p^ ) = \/x. By contrast to these 
simple formulas, there is no analytical formula for the REE 
for p( A > for arbitrary a and p, except some special cases. For 
example, the REE can be calculated analytically, according to 
Eq. ( fTTT ). for the Horodecki states, given by Eq. (|26V 

Our numerical calculations shown in Fig.[TJc) indicate that 
the Bell-diagonal states are likely to be the lower bound for 
the REE vs nonlocality. Thus, let us analyze whether the ADC 
states can be diagonal in the Bell basis. By denoting the Bell 
states as follows |ft) = \ip~), |ft) = \ip + ), |ft) = \<f>~), and 
|ft) = |</> + ), we find that the ADC state can be given in the 
Bell basis as: 

p^ (a, ft, g 2 ) = r_|ft > (ft | + r + |ft) (ft | + r|ft) (ft I 

+r|ft) (ft| + («|ft) (ft | + r|ft) (ft | + h.c), (29) 

where t = a(l - q 2 ) + r - \ ; r± = i[l - 2r ± C(p^)}, 
and r = |[(1 — a)q\ + aq 2 ]- So, this state is diagonal in the 
Bell basis only if a = and 1 . Note that for a = 1/2, one has 
r_ = \{pi — p 2 ) 2 , which vanishes for p\ — p 2 . but another 
off-diagonal term r — \{q\ + q 2 ) vanishes only if there is no 
damping at all. This shows (as also confirmed numerically in 
Fig. |2) that the ADC states do not have the minimum of the 
REE for a given nonlocality except only two points for B = 
and 1. By contrast, Figs.[Tfc) and|2]show that the ADC states 
are likely to be the upper bound of the REE vs nonlocality. 
This observation is confirmed analytically in Sec. [V] 



B. Phase-damping channel 

The Kraus operators for the single-qubit phase-damping 
channel (PDC) can be given by y[]: 

#ofe) = |0)(0| + ^|l)(l|, E 1 (q i ) = V«11><1|, (30) 

where qi = 1 — pi is the phase-damping coefficient. Thus, 
one can find that the pure state \t/j a ) after being transmitted 
through the PDC is changed into the following mixed state 

p^ dc \a, gi ,q 2 ) = o|01><01| + (1 - a)|10)(10| 

+y(|01)(10| + |10)(01|), (31) 

where y = -J a(l — a)p\p 2 . This state can be also given as 

p^ dc \a, qi , q2 )=p"\M(M+q"\0l)(0ll (32) 

where the effective damping constant can be defined as q" = 
1 — p" = a(l — pip 2 ) and the pure state \tp a //) is given by 
Eq. (l2Tb but for a" — apip 2 /p" ', which takes the values < 
a" < a. Note that if a = Pip 2 /(qi +q 2 —qiq 2 ) then a" = a. 

We find that the nonlocality, negativity and concurrence for 
the PDC state are equal to each other and are given by 

B( p (pdc)) = jV( p (pdc)) = C( p (pdc)) = 2 y, (33) 

The range of the REE for a given nonlocality for the PDC 
state, given by Eq. Pit , is shown in Fig. [2] Note that it is 
unlikely that there is analytical formula for the REE for the 
PDC states for arbitrary parameters a, q\, and q 2 . However, 
the REE can be found in some special cases. For example, if 
a = [1 + = a 1 then Eq. ( l32b reduces to 

p< pdc V>9i.«») = P {V) = 2(1 - a')|^ + )(V + l 

+(2a'-l)|01)(01|, (34) 

which means that becomes the Bell state \ip + ). For this 
state, one finds that 

B' = B(pW) = C{pW) = A^(p (V) ) = 2(1 - a'). (35) 

The state p^\ the same as pure states, reaches the upper 
bound for the negativity versus concurrence lli5ll . The CSS for 
reads as a( v ) = (1 - B'/2)|01)(01| + (B'/2)|10)(10|. 
Thus, the REE can simply be given by: 

E R {p^)) = h{\B>) h{\W{l-B'f + {B'f + l]}. 

(36) 

Now, we show that the minimum of the REE vs nonlocality 
can be reached by the PDC states. In the Bell basis, the PDC 
state is given by 

p^\a, qil q 2 ) = (i - y)|ft)<ft| + (§ + y)|ft)(ft| 

+(a-|)(|ft)(ft| + |ft)(ft|), (37) 

which clearly becomes diagonal for initial Bell states (i.e., for 
a = 1/2), given by Eq. (fT8l for q = l—p= 1/2 — y (see e.g. 
Ref. Jll). The REE for the state p (D) is given by Eq. (|20b . 
Hereafter, we consider only the case, when a — 1/2 and use 
the shorthand notation p( D ) = p( D2 ) . 
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V. KARUSH-KUHN-TUCKER CONDITIONS FOR THE 
REE VS NONLOCALITY 

In this section, we derive the KKT conditions, as a general- 
ization of the method of Lagrange multipliers, in order to find 
the states with extremal values of the REE for a fixed value of 
the nonlocality. 

Thus, let us consider the following Lagrange function: 

C = Tr(plogp) - Tr(p log a) + l[Tr(pBcHSn) - (3} 

- Ti(Xp) + A(Trp — 1), (38) 

where I, A, and X are Lagrange multipliers and (3 — 
2yM(p). For a small deviation of p — > p + A, we have 

C -> C + Tr(A log p) - Tr(A log a) 
1 1 



from which I is determined. Then, there must exist X such 
that 



Trip 



-A- 



-dz 



lo P + z p + z 
+ ZTr(AScHSH) - Tr(AX) + ATrA 
C + Tr[A(log p - log a + P 

+ ZBchsh-X + A)], 



(39) 



where P is a projector on the support space of p. Thus, the 
KKT conditions on the parameters I, A, and X are the follow- 
ing: 

logp-lo g( 7 + P + /Z?cHSH-* + A = 0, (40a) 
X > 0, Tr(Xp) = 0. (40b) 

From the condition (140al>, we obtain 



(41) 



= Tr(plogp) - TrOologcr) + Tr(pP) 
+ /Tr(p6 CH SH) - Ti(pX) + XTip 

= jE R (p) + l + l(3 + \, 



where 7=1/ log 2 e and e is the base of the natural logarithm. 
Thus, A = —-fEn(p) — 1 — 1/3. Now, we can rewrite Eq. (140a) 
and obtain the simplified KKT conditions: 

= logp - logo- + P + ZSchsh - X 

-lE R {p) -1-W, (42a) 
X > 0, Tr(Xp) = 0. (42b) 

One can easily confirm that the condition 

P [ log p - log a + ZBchsh - 7^i? (P) - W] P = (43) 

is a new KKT condition since FA = 0. 

The boundary states are rank deficient. We search for them 
among the rank-2 states. We have a strong numerical evidence 
(see points for randomly-generated density matrices in Fig.[T) 
implying that the extreme two-qubit states can be found within 
the rank-2 class of mixed states. Thus, let us now consider the 
case where p is a rank-2 mixed state, i.e., p = Ai|ei)(ei| + 
A2|e2)(e2|, where A^ are non-zero eigenvalues of p, and | e^) 
are the corresponding eigenstates. From Eq. (142b) . we have 

(ei| log<7|e 2 ) = /(e 1 |ScHSH|e 2 >, (44a) 
jE R (p) +1(3 = log Ai - (ei I log cr|ei) 

+ /<ei|B C HSH|ei), (44b) 



X = log p - log a + P + IBchsk 
-[ 7 E R (p) + l + lp] >0, 
Tr(Xp) = 0. 



(45a) 
(45b) 



The KKT conditions, given by Eqs. (44ah — f45bl) . are the nec- 
essary conditions in searching for the boundary (extreme) 
states among rank-2 mixed states. 



A. Lower bound of REE vs nonlocality 

Here, we show the KKT conditions are satisfied by the Bell- 
diagonal states, given by Eq. ( fT8l , as obtained by the phase 
damping of the Bell states \ip + ) or, equivalently, assum- 
ing pf> 1 — p. 

For these states the CSS is given 

^ (D) = ^ + )(^ + l + lVO<V^I), (46) 
and the Bell-CHSH operator is 

S chsh = v[ - CT 3 ® CT 3 + (2p - 1)<ti ® en] 

= 27 ? [p(i^ + )(^+i-i<r)(<ri) 

+ (l-p)(|^-)^-|-|^+)^+|)], (47) 



where -q = 2/y/l + (2p — l) 2 . So, by applying Eqs. © 
and (O, one can find the simple expression, given by Eq. ( fT9l , 
for the nonlocality. 

Since it holds (ei| logcr( D )|e 2 ) = (ei|g^, H [ea) = 0, 
Eq. (144a) is satisfied for any I. From Eq. ( 144b) . / is deter- 
mined through 



logp + 1 + 2rjpl = 7 ^ D) + ip, 



where (3 — 2 yT+ (2p — l) 2 and the REE is given by 

Eg ) =E R (pW) = l-h(p) ) 



(48) 



(49) 



where h(p) is the binary entropy defined below Eq. (|8)- Then, 
the KKT conditions are satisfied if 

X = - [ 7 i4 D) + 1 + ip + 2r,l(l - p)]|0+> (0+ 1 

- [ 7 i4 D) + 1 + 1(3 + 2rjlp] \<j>- ) {(f)- 1 > (50) 

for all values of p. Hence, the Bell-diagonal states can yield 
the extreme (in this case minimum) value of the REE for a 
given value of the nonlocality. This conclusion is confirmed 
by our Monte Carlo simulations (see Fig . [TJ . 



B. Upper bound of REE vs nonlocality 

Here, we show that the KKT conditions are also satisfied 
by the amplitude-damped states p' A ' (a,p), given by Eq. d25l l, 
for properly chosen parameters a and p. These MEMS having 
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the highest REE for a given negativity are shown by curves A" 
inFigs.[Hc)and[lb). 

The support space of p^ (a, p) is given by the eigenvectors 
|ei) = 1 00) and |e 2 ) = \ijj a ). The corresponding CSS is given 
as Hi: 



a {A) (a,p) = i?i|00)(00| + i? 4 |H)(H| 

+A + |A + }<A+|+A_|A_><A_|, 



(51) 



where 



A+ = x 



Ri 



i?3 ± \/ (i?2 — -R3) 



4i?ii? 4 



|A±) = A± [(A± - R 3 )\Q1) + y/RilhllO) 

normalized by A± = [(A± - i? 3 ) 2 + RiR 4 ] ~ 1/2 
late the CSS we have to compute 



(52) 
To calcu- 



R 2 = 


1 
4 


1 + 3(1 -p) + 2pa -ARi-VS 


1 


i?4 = 


Ri 




(53a) 


R 3 = 


1 - 










i=l,2,4 




5 = 


(4 


- ipf - 4a(l - a)p 2 - 8i?i(2 - 


P) 



+16y / Ri(R 1 - l+p)p 2 a(l - a), 
where Ri is obtained by solving 

ap = R 2 + 2R 4 (Rl - R 2 R 3 + 2R 1 R 4 )/z 2 



-2R 4 {R 2 -R 3 )/(Lz) 



(53b) 



(54) 



where z = ^{R 2 - R3) 2 + AR 1 R 4 and L = log(i? 2 + R3- 
z) — log(i?2 + i?3 + z). The Bell-CHSH operator in this case 
reads as 



g(A) 
°CHSH 



where 771 



171 



■\2 



(1 - 2p)a® 2 + 2p^{l - a)aaf 
if4p 2 (l-a)a- (l-2p) 2 < 0, 



V2 p^/(l-a)a(af 2 +af 2 ) 
otherwise, 



(55) 



2/^(1 - 2p) 2 + 4p 2 a(l - a), m 



2/\/8p 2 (l — a) a. Thus, by applying Eqs. © and (|6), we 
find that the nonlocality is given by Eq. d27b . Since we fix the 
value of p or equivalently B = B(p( A >), we can express p in 
terms of B and parameter a in the following way 



'l(2p- VWP 2 -B 2 -1) 



if 2V2 + 2S 2 < 4p < 2 + V2 + 2B 2 , 
^(P~V5P 2 ~4P-B 2 ) 

if 2 + V2 + 2B 2 <Ap<A. 



(56) 



We can easily check by using Eqs. d55l ) and ( BTT ) that the con- 
dition, given by Eq. d44ab . is always satisfied. Thus, the condi- 
tion, given by Eq. d44bb . is also satisfied. To check the remain- 
ing conditions, we need an explicit expression for the REE 
and multiplier I, which is equivalent to solving Eq. (|54l . This 
equation contains logarithms and can be easily solved only in 



special cases as, e.g., the Horodecki states (a — 1/2). Never- 
theless, our numerical analysis reveals that for p( A \a,p), for 
which the REE reaches the largest value for a given value of 
nonlocality, Eq. (I45bb is satisfied for the following coefficients 
of the amplitude-damped states: 



\ (2 + V2 + 2B 2 ) if B<B Q , 
1 HB>B , 



(57a) 



1 



1>P> 1 (2+ \/2 + 2S 2 ) ifB=B , (57b) 



a = — ( p — v hp 1 — Ap — B 2 
2p \ 



(57c) 



where B = 0.81686 [see Fig.Etb)]. 

Thus, by comparing the REEs for a given nonlocality for 
the optimal mixed states (denoted by subscript A") and pure 
states, we can conclude that 



E {A " ] {B) > E R P) {B) for < B < B 



J R 

,(A"), 



(58) 



E£ '(B) < E^'(B) for B Q < B < 1, 



where B Q = B{p 4 ) = B(p 5 ) = 0.81686 and E (A "\b ) = 

E [P) (B ) = 0.7445 [see Tab. I and Fig. 3(b)]. On the other 
hand, by comparing the REE for a given negativity for the op- 
timal mixed states (denoted by subscript A') and pure states, 
it holds Ql: 

E R A '\N) > E R P) (N) for < N < N , (59) 
E (A '\n) < E R P) (N) for JV < AT < 1, 

whereto = N(p 2 ) = 0.5271 and E R A ' ) (N ) = E { r \n ) = 
0.3847 [see Tab. I and Fig. Oa)]. As a reminder, subscript 
A' (A") corresponds to the amplitude-damped states p^ A ) 
(p( A )) optimized for the REE vs negativity (nonlocality). 
As shown in Fig. [T] the states p( A ' and p( A > are, in gen- 
eral, inequivalent. The maximal differences, as defined by 
Eq. QUI, are equal to AE R (N) = 0.2018 for N = 0.7071, 
and AEr(B) = 0.4040 for B = 0. It is seen that ranges of 
mixed states, which are more entangled than pure states, are 
much smaller for the REE vs negativity in comparison to the 
REE vs nonlocality. 



VI. CONCLUSIONS 

We studied the relation between the relative entropy of en- 
tanglement Er and the nonlocality measure B correspond- 
ing to a degree of the violation of the Bell-CHSH inequality 
in two-qubit systems. We found states of the extremal value 
of Er for a given value of B. We showed that the obtained 
states satisfy the Karush-Kuhn-Tucker conditions, derived in 
Sec.|VJ as well as they provide a boundary for the REE values 
obtained by our Monte Carlo simulations presented in Fig. [TJ 

We demonstrated that mixed states can be more entan- 
gled in terms of En than pure states if the nonlocality B E 
(0,0.82) and the maximal difference between these REEs is 
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AEr = 0.4 as shown in Fig. [2b). As discussed in Ref. lfl7ll . 
Er as a function of the negativity N can also exhibit this prop- 
erty but (i) the maximal difference AEr is one-order smaller 
(i.e., AEr = 0.039) and (ii) mixed states are optimal for a 
shorter range of the negativity, i.e. N 6 (0, 0.53), as pre- 
sented in Fig. [3f a). For appropriate comparison, we normal- 
ized B to be equal to N for any two-qubit pure state. 

We showed that these maximally-entangled mixed states 
can be obtained from pure two-qubit entangled states by lo- 
cally subjecting one or both qubits (of the entangled pair) to 
the amplitude-damping channel either. While the minimally- 
entangled states can be generated from pure states by sub- 
jecting one or both qubits to the phase-damping channel. We 
found that the amplitude-damped states (yellow and green ar- 
eas in Fig. [2]) are more entangled than the phase-damped states 
(blue and green areas) for a given nonlocality B < 0.8169 
(corresponding to the nonlocality of the states p4 and p$ 
shown in Fig.Q]and Tab. [I]). 

However, for values B > 0.5856 (point pe) there exists 
a range of states, obtained either by the amplitude or phase 
damping, which have the same value of En for a given B 
(green area in Fig. |2). Moreover, we found that the upper 
bound on the REE of the phase-damped states (cyan and green 
areas in Fig. |2|i for a given nonlocality is provided by pure 
states (curve P in Fig.|2|. However, pure states have the high- 



est REE only for B > 0.8169. 

Thus, we found that for a large range of the nonlocality, 
mixed states can be more entangled that pure states, and, sur- 
prisingly, these mixed states can by obtained by the ordinary 
amplitude damping of pure states. 

We believe that these results can stimulate further quest 
for practical protocols of quantum information processing for 
which mixed states are more effective than pure states. 
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